Zeeman ratchets for ballistic spin currents 
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We investigate the possiblity of creating directed spin-polarized currents in a two-dimensional 
electron gas (2DEG) subject to an asymmetric magnetic field and an external adiabatic driving. We 
thereby generalize concepts of quantum charge ratchets to the case with spin. Due to the Zeeman 
term in the Hamiltonian, spin-up and -down electrons experience different effective potentials, which 
can be tailored to achieve net spin currents without corresponding charge currents. We consider 
ballistic, coherent transport in waveguides defined in a 2DEG, where the magnetic field modulation 
is induced by ferromagnetic stripes on top of the 2DEG. 
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I. INTRODUCTION 

Magneto- and spin-electronics ■ represent prominent 
fields of modern condensed matter physics. The vi- 
sion is to create electrical and opto-electronic de- 
vices with new functionalities based on spirA Since 
semiconductor-based spintronics requires the generation 
of spin-polarized currents, various proposals have been 
made to this end, ranging from spin injection from a fer- 
romagnet to adiabatic spin pumping 3 . 

Here we propose an alternative way to generate spin 
currents, namely by generalizing ratchet mechanisms for 
creating directed particle currents^ to the case of spin. 
Particle ratchets are based on spatially periodic po- 
tentials with broken inversion symmetry which - upon 
ac-driving with zero net bias - give rise to a directed 
particle flow. They have been realized in ballistic mi- 
crostructures^. In the present Paper we suggest to tailor 
spatially-periodic, non-uniform magnetic fields such that 
spin-up and spin-down particles exhibit opposite directed 
motions due to the corresponding Zeeman potential, re- 
sulting in a net spin current. A corresponding proposal 
for a spin ratchet based on spin-orbit interaction has been 
recently worked oulpi. 



II. MODEL AND FORMALISM 

We consider a two-dimensional electron gas (2DEG) 
in the (x, y) plane and a quantum wire in the a>direction 
connected to two non magnetic leads. The system is em- 
bedded in a non-uniform magnetic field B(x) — B y (x)y 
pointing in the y-direction. The Hamiltonian of the sys- 
tem reads 
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B y (x)a y + V(y), 
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where g* is the effective gyroscopic factor, m* the effec- 
tive mass, /ib the Bohr magneton and a v the Pauli spin 
operator. V(y) denotes the lateral confining potential. 

The particular configuration of the magnetic field per- 
mits to neglect orbital effects and allows for solving the 



time-independent Schrodinger equation through the sep- 
aration ansatz: $ CT (x,j/) = <p a {x)x n {y)W)y Here Xn(y) 
are the transversal eigenfunctions, <p a (x) are plane waves 
and \a) y are spin eigenvectors along the y-direction. The 
specified magnetic field B(x) can be realized by making 
use of the fringe fields of ferromagnetic stripes patterned 
on top of the 2DEG and magnetized in the in-plane di- 
rection perpendicular to the quantum wire. The mag- 
netic field of a stripe of length b along the y-direction 
(with b greater than the other two spatial dimensions) 
is in its central region (y w 0), almost homogeneous in 
y-direction, B y 3> B x , B z . Then the expression for the 
magnetic field in Eq. is valid if the quantum wire has 
a width W <C b and lies right beneath the center of the 
stripe (y rj 0). In this limit the other components, B x 
and B z , are small compared to B y and their effect to mix 
eigenstates \a) y of the system can be neglected. 

We are interested in the quantum transport properties 
of this system within the Landauer approach. For the 
Hamiltonian the spin-flip transmission and reflection 
probabilities are T CTi _ ff = T' a _ a = R a - a = R' a .-a = 
due to the suppressed mixing of the spin channels^. Fur- 
thermore, different lateral modes are decoupled and as 
consequence it is sufficient to look at transmission and 
reflection probabilities of the lowest energy sub-band 
(n = 1) only, since the m-th mode will just mimic this 
behavior with an energy offset of E m — E\ . 

In order to induce particle motion between the two 
leads, we consider an adiabatic square-wave driving as 
a rocking potential. The system is periodically switched 
between two rocking conditions, labeled by ±f/o {Uq > 
0), where it can be considered in a steady state. There- 
fore, the averaged net charge and spin currents inside the 
wire read 



(ic) = ~[W+u ) + W-u )] 



and 



-[IU+Uo)+IU-Uo)} 



(2a) 



(2b) 



respectively. 
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FIG. 1: (left panel) Spin-up T+,+ (dashed line) and spin-down Tl,- (dotted line) transmissions and their difference (solid 
line) as a function of the injection energy for a magnetic profile (|1J with N = 1, a — 0.22 and Uo = 0.09 Ub in a fixed rocking 
situation +Uo- The curves for the opposite rocking situation — Uo can be obtained through the relation (jjj. The Inset shows 
the modification of the effective Zeeman potential barrier due to the voltage drop in the two rocking situations. 

FIG. 2: (right panel) Spin-up (dashed lines) and spin-down (solid lines) transmissions as a function of the injection energy in 
the two rocking situations for the magnetic profile Q with N = 1.5, a = 0.22 and Uo = 0.1517b. The upper (lower) Inset 
shows the difference of spin-up (spin-down) transmission in the two rocking situations. 



In steady state these currents are obtained within the Landauer formalism as: 

J c = -^/ dE[f(E;^)-f(E;^ R )][T +!+ (E)+T_ t _(E)] , (3a) 

i r°° 

% = ^J o dE [/(^ Ml) - f{E; m )} [T +I+ (E) - T--(E)] . (3b) 



In Eqs. ©, T + ^ + (E) and T--(E) are the spin resolved 
transmissions, and /(.E; /Hl/r) is the Fermi function for 
the left /right lead at chemical potential //l/r- The ab- 
sence of processes that can change the spin state of the 
electron causes the spin current to be conserved through- 
out the wire. 

Charge and spin currents © do not vanish in the non- 
linear regime. To model a finite voltage drop across the 
two leads, we add the term Hu(t) = —U(t)g(x,y; U) to 
the Hamiltonian (JTJ. Here the function g(x,y;U) de- 
scribes the spatial distribution of the electrostatic po- 
tential inside the mesoscopic system and in general is 
obtained through a self-consistent solution of the many- 
particle Schrodinger equation and the Poisson equation. 

In the present Paper we employ a heuristic model for 
the voltage drop g(x, y; U)^. We assume that in the cen- 
tral part of the mesoscopic system the voltage drops lo- 
cally proportional to the gradient of the magnetic field 
■§^g(x, y\ U) oc | -g^B y {x)\. This model is justified by the 
fact that the Zeeman term acts as an effective potential 
barrier, and it has been shown£ that a more rapid poten- 
tial variation leads to stronger wave reflection and hence 



to a more rapid local voltage drop. As a direct conse- 
quence, the motion of spin-up and spin-down electrons 
can give rise to different rectified currents, opening the 
possibility for the spin ratchet effect. 



III. NUMERICAL RESULTS 

The energy dependent transmissions T+.+ and T_,_ 
are numerically calculated from lattice Green functions. 
Within the scattering region of length L and width W, an 
array of magnetic stripes with antiparallcl magnetization 
can be chosen such that the resulting magnetic field is 
given by 



B y (x) = B Q 



— sin 



2tt \ /4tt 
— x j + a sin I —x 



(4) 



Here I is the spatial period of the non-uniform magnetic 
field, and a is an asymmetry parameter. The ratio N — 
L/l is the number of magnetic field unit cells within the 
scatterer. The system parameters are chosen to have 
the experimentally accessible values^: m* = 0.1m, g* — 
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FIG. 3: (left panel) Averaged spin-up (solid line) and spin-down (dashed line) current as a function of the Fermi energy for a 
magnetic profile with N = 7.5, a = 0.22 and U = 0.09 Ub- 

FIG. 4: (right panel) Averaged spin (solid line) and charge (dashed line) current as a function of the Fermi energy for a magnetic 
profile with N = 7.5, a = 0.22 and U = 0.09 U s . 



40, W = lOOnm and I = 2/Ltm. The magnitude of the 
magnetic field Bq is fixed to 0.1T. 

The system Hamiltonian {Q, in presence of the mag- 
netic field Q and A~ integer, shows the symmetry rela- 
tions [H,a x R x ] = and [H,a z R x ] = 0, where R x is the 
operator for reflection at the x-axis. Those allow us to 
interrelate the elements of the scattering matrix for the 
two rocking situations. We find a non-zero averaged spin 
current 

TaA+ u o) = T.^i-Uo) => (7|") = m+U ) (5) 

and vanishing averaged charge current (7c) = 0. 

In Fig. ^ we present for the case of N = 1 the trans- 
mission probabilities as a function of the injection energy 
normalized to the maximum height of the magnetic bar- 
rier, Ub — 9*(J-b xasx[B y (x)]/2. For a fixed rocking situa- 
tion, spin- up and spin-down electrons pass through oppo- 
site magnetic barriers with a different maximum height 
due to the effective voltage drop (see Inset of Fig. [Q. 
This causes the difference between the T +i+ and T_._ 
transmissions. Apart from this classical effect, this trans- 
mission difference can also have a quantum origin, e.g. 
tunneling, relying on the asymmetric shape of the barri- 
ers (a ^ 0). To exclude the classical effect, we turn to a 
magnetic field profile where the same voltage drops from 
the high potential reservoir to the first barrier in either 
of the two rocking situations. This can be achieved by 
extending the size of the scattering region in by half 
a unit cell, i.e. N = 1.5. 

Since N is no longer an integer, Eq. JSJ does not hold 
true anymore. Rectification of spin-up and spin-down 
electrons still may provide us with non-zero averaged 
spin current, but due to the absence of symmetry the 
averaged charge current does not vanish anymore. In 
Fig. [21 we show the energy dependence of the transmis- 



sion probabilities for both spin species in each rocking 
situation and in the Insets their difference for A" = 1.5 
and a = 0.22. The main difference between T + + (black 
dashed lines) and 7 1 -,- (red solid lines) is due to the pecu- 
liar choice of the magnetic field profile: spin-up electrons 
travel through a single barrier whereas spin-down elec- 
trons experience a double barrier structure. Indeed, this 
double barrier structure gives rise to resonant tunneling 
behavior for T_ _. Due to this fundamental difference 
it is intuitively clear that the rectification for both spin 
species is not opposite to each other and therefore results 
in finite averaged charge and spin currents. 

Increasing the number of unit cells of the magnetic 
profile and thereby approaching the setting of a periodic 
ratchet, it is possible to observe how the resonances are 
giving rise to the formation of minibands. In Fig. [3] wc 
show the averaged spin resolved currents as a function of 
the Fermi energy ^o- Those are evaluated for a magnetic 
profile with A" = 7.5, a = 0.22 and the magnitude of the 
applied voltage is Uq — 0.09 Ub- It is evident that for 
almost all the values of in the plot the averaged spin- 
up current has opposite direction and approximatively 
equal magnitude with respect to the spin-down current. 

In Fig. 0] we show the averaged spin and charge cur- 
rent as a function of the Fermi energy fj,Q for the same 
parameters as in Fig. |3J The averaged charge current 
(dashed line) is mostly small compared to the averaged 
spin current (solid line). The sharp peaks and the oscil- 
lating structures in Figs. El and 0] can be attributed to the 
minibands of the magnetic superlattices. 

Those results confirm our expectation that the asym- 
metry in the barriers of the magnetic field can translate 
into different, sometimes even opposite rectification for 
spin- up and spin-down electrons. 
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